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Abstract 

We consider a network model where the nodes are grouped into a number of clusters and propose a 
distributed dynamic frequency allocation algorithm that achieves performance close to that of a central- 
ized optimal algorithm. Each cluster chooses its transmission frequency band based on its knowledge 
of the interference that it experiences. The convergence of the proposed distributed algorithm to a 
sub-optimal frequency allocation pattern is proved. For some specific cases of spatial distributions 
of the clusters in the network, asymptotic bounds on the performance of the algorithm are derived 
and comparisons to the performance of optimal centralized solutions are made. These analytic results 
and additional simulation studies verify performance close to that of an optimum centralized frequency 
allocation algorithm. It is demonstrated that the algorithm achieves about 90% of the Shannon capacities 
corresponding to the optimum/near-optimum centraUzed frequency band assignments. Furthermore, we 
consider the scenario where each cluster can be in active or inactive mode according to a two-state 
Markov model. We derive conditions to guarantee finite steady state variance for the output of the 
algorithm using stochastic analysis. Further simulation studies confirm the results of stochastic modeUng 
and the performance of the algorithm in the time-varying setup. 

I. Introduction 

Dynamic frequency allocation has an important role in the performance of wireless ad-hoc 
networks, for it results in less transmission power, which is a crucial objective in communication 
networks. To do this in an optimal way, there needs to be a centralized processor with full 
knowledge of the spatial distribution profile of the network clusters. However, in many emerging 
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wireless networks (such as ad hoc networks, cognitive radios, etc.), no central frequency allo- 
cation authority is naturally available. This makes distributed frequency allocation an important, 
but mostly unchartered territory in wireless networking. 

Centralized frequency allocation has been extensively studied in the literature (Please see [8] 
and [11]). There are a number of proposed solutions to similar problems in different contexts 
(Please see [2], [4], [6], [9], [10], and [12]). These include methods based on graph coloring 
for cognitive networks, iterative waterfilling for Digital Subscriber Lines (DSL), etc. These 
approaches either excessively simplify the interference models, or are not fully decentralized, or 
require too much information exchange between autonomous nodes/clusters, or suffer from all 
these shortcomings. Additionally, they are all too complex to implement. In [10], the approach is 
based on approximating the optimal resource allocation solutions on a graph. Others [9] propose 
that secondary users choose their spectrum according to their information about their local 
primary and secondary neighbors. They employ a simplified model for mutual interference of the 
network nodes that turns the problem into the graph multi-coloring problem. They subsequently 
compute a sub-optimal solution to the graph multi-coloring by using an approximation algorithm 
to the graph labeling problem. 

In the context of Digital Subscriber Lines (DSL), some recent works regarding spectrum 
balancing have been done (Please see [2] and [12]). The objective of spectrum balancing in DSL 
systems is to maximize the throughput of each user by shaping its Power Spectral Density (PSD) 
of transmission, satisfying a certain power constraint. In [12], a method of iterative waterfilling 
is proposed in order to solve the problem. In the case of two users, they show the existence and 
conditions on the uniqueness of a Nash equilibrium point for the iterative algorithm. However, 
each user must know a weighted sum of the PSD of the other users (interference), in order to 
do waterfilling. The iterative waterfilling algorithm has high complexity and the resulting Nash 
equilibrium point is not necessarily the optimal solution. For instance, in a two-user scenario, if 
both users start with a flat PSD initially, iterative waterfilling does not change their PSD. This 
is clearly a Nash equilibrium point, but is far away from the optimal answer. In [4], it is shown 
that this non-optimal Nash equilibrium point might be the only Nash equilibrium, and therefore 
iterative waterfilling fails for various scenarios. 

In [2], the users need to balance their power along a number of tones in order to optimize their 
throughout under power constraints. The optimization problem is relaxed based on introducing 
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a virtual user with fixed thresholds. It turns the problem into a separable optimization problem 
across the tones for different users. An algorithm has been proposed to solve the relaxed 
problem iteratively via solving local optimization problems by the users. The knowledge of 
a weighted sum of the PSD of the other users (interference) is required for each user to solve 
its local optimization problem. The convergence of the algorithm has been shown in high SNR 
regime. Simulations show that the achievable region resulted by the solution of the relaxed 
distributed optimization is close to that of the optimal centralized solution. However, no one-to- 
one correspondence between the points of the achievable regions of the optimal (centralized) and 
decentralized algorithms is guaranteed. Therefore, the algorithm does not necessarily converge 
to optimal values. For the case of asynchronous transmission (in the presence of ICI), the 
optimization problem is not separable across the tones. They have therefore used heuristic 
optimization approaches with no convergence guarantees. 

In [4], it is shown that the problem of optimal PSD shaping across the users is reducible 
to that of allocating piece-wise constant powers. This result reduces the complexity of the 
spectrum sharing problem. Furthermore, a number of achievability and existence results in the 
context of non-cooperative and cooperative game theory for obtaining efficiency and fairness have 
been established. Another approach has been presented in [6], where each user in the network 
announces a price to the other users to adapt the power allocation accordingly. Convergence 
results have been established using supermodular game theory. 

In the model we consider, the nodes are divided into different clusters and each cluster is 
represented with a cluster head. This is motivated by the fact that various networks are naturally 
clustered (e.g., combat scenario, WLAN Hotspots, WPAN). Each cluster head, having knowledge 
about the interference it experiences, chooses the frequency band with the least amount of 
interference from the other clusters. The channel model we consider is the common path loss 
model, which gives a more refined model than that of the existing literature. It is shown that 
this distributed strategy converges to a sub-optimal spectrum assignment, without any cross- 
cluster information exchange. In other words, the algorithm converges to a local minimum of the 
aggregate interference of the network. Simulation results (Section jV]) show that the minimization 
of the aggregate interference of the network results in a sub-optimal solution for the problem of 
maximizing the aggregate Shannon capacity of the network links. It must be noted that regardless 
of the model, channel reciprocity is sufficient for the convergence of the algorithm. 
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The proposed algorithm provides a simple, fully distributed, dynamic frequency allocation 
strategy that requires neither any information exchange between autonomous devices, nor even 
any knowledge of the existence of other autonomous entities. Additionally, the proposed algo- 
rithm can be used in conjunction with any realistic wireless radio channel model such as those 
commonly employed in wireless standards (Hata model, Okumura model, etc.). We will further 
propose performance bounds on this sub-optimal spectrum assignment for some specific network 
topologies. 

We also present a framework to analyze the scenario that clusters can be in sleep or active 
mode, and go off and on according to time- varying statistics. In our model, the activity of the 
clusters is described by a stochastic process. For the simplicity of the analysis, we assume that all 
the clusters go on and off according to a two-state Markov model independently. We have shown 
that the temporal dynamics of the algorithm can be accurately modeled by an exponential decay. 
We have derived a stochastic dynamical equation for the evolution of the aggregate interference 
of the network. The time variation of the activity of clusters is also included in the model, which 
results in a stochastic differential equation for the steady state behavior of the algorithm. We 
have further derived a trade-off inequality in terms of the update rate, switching rate between 
sleep and active mode and the geometrical properties of the distribution of users, guaranteeing 
a certain steady state variance on the performance of the algorithm. Simulation results verify 
the accuracy of stochastic modeling and the performance of the algorithm in the presence of 
time-varying statistics. 

The outline of this paper follows next. In Section UIl we discuss the system model and the 
assumptions employed for our analysis. We disclose the Main Algorithm in Section|IlIl In Section 
HVl we present the main results of this paper regarding the convergence and performance of the 
main algorithm. Simulation results for some specific network topologies are provided in Section 
IVl Finally, we discuss the contributions of the proposed method and future areas of research in 
Section |VIl The proofs of main results can be found in Appendices IJ and |Ill 

II. System Model 

Suppose that we have a set of network nodes distributed in space such that they can be 
partitioned into a union of clusters. These networks often happen in nature, for instance, in a 
combat scenario a group of soldiers can be divided into a number of clusters according to their 
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missions. Communication within clusters is then very desired. There are a number of efficient 
methods for partitioning the network elements, but this topic is not the focus of this paper. We 
assume that the clusters are already formed in a specified manner. 

In light of the above, our network model is given by collection of nodes in clusters, q, 
i = 1, - ■ ■ , N, where each cluster has a cluster head responsible for managing some of the 
network functions. Let dij denote the distance between the cluster heads of q and Cj (Fig. 1). 
We use the following assumptions for the model: 

• The ith cluster, q, contains Ui users. 

• At each time slot for any cluster at most one user is transmitting and one user is receiving. 
For simplicity of mathematical analysis, we have assumed that at each time slot for any cluster 
exactly one user is transmitting and one user is receiving (See Sections IIV-FI and jV]). This 
assumption can be relaxed to any scenario satisfying channel reciprocity between clusters. For 
instance, in an alternative scenario each user transmits and receives its data through the cluster 
head. This model also satisfies the channel reciprocity conditions. Therefore, the results can be 
generalized to various other models. 

• Each user transmits with power KPq, where is a constant such that the power at 1 meter 
is Pq. The assumption of equal transmission powers can be further relaxed and is adopted for 
mathematical convenience. 

• The distances between clusters are much larger than the size of clusters and bounded below 
by a distance 6. 

• The transmission model is path loss with exponent i]. No shadowing and fading is assumed 
in the current analysis. However, the analysis can be generalized to more realistic models of 
transmission. 

• The accessible spectrum is divided into r different bands, denoted by 6i, ■ ■ ■ , &r- 

• At time t, the ith cluster is in state Si{t) G {1, 2, ■ ■ ■ , r}, corresponding to the index of the 
transmission band it is using. 

• The probability of two clusters updating their frequency bands at the same instance of time 
is negligible. This assumption can be relaxed [1]. 

• The rate of change of the spatial distributions of the clusters in the network is much less 
than the processing/transmission rate. Therefore, the topology of the network is assumed to be 
fixed in the analysis of the frequency allocation algorithm. 
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• We choose the aggregate interference of the network as our performance metric. This 
assumptions makes the analysis mathematically tractable and is a reasonable metric for perfor- 
mance comparison. Simulation results in Section |V] show that the minimization of the aggregate 
interference of the network results in a near-optimal solution for the problem of maximizing the 
aggregate Shannon capacity of the network links. 

III. The Algorithm 

Using the above assumptions, we approximate the ith cluster, q, by a single node with 
transmission power Pq within a distance dij from the other clusters. The interference experienced 
by Ci caused by all the other clusters is therefore, 

I,^{N,{d,,},t) = J2^6is,it),s,it)) (1) 

where 5{x,y) is the Kronecker delta function, defined as 

1 X = y 

S{x,y) = { (2) 
[0 x^y 

At time t = tn, one of the clusters, say q, updates its transmission frequency band. The 
nature of the update procedure is asynchronous for all the clusters. This is intuitively appealing, 
because of the nature of ad-hoc networks, where there is usually no common clock among the 
nodes. We assume that the updates are taking place at times ti,t2,---. We can therefore change 
the continuous-time interference model in Eq. ([T]) to a discrete-time version as 

4(iV,{d,,},/) = 5^^5(s,(/),s,(/)) (3) 

where / corresponds to the time t = ti, when an update is taking place. Let Ck{N, {dij}, I) denote 
the set of clusters transmitting in band bk prior to time /. Also, let I^^^kiN, {dij},l) denote the 
interference experienced by Cj caused by all the clusters in Ck{N, {rfy }, I) if Cj was transmitting 
in band bk, for k = 0, ■ ■ ■ ,r — 1. Ic^^kiN, {dij}, I) can be written as 

/,,,,(iV,{d,,},/)= J2 
We denote the aggregate interference of the network at time / by I {N, {dij}, I) as 

i{N, {d.,}, i) = j2 UN, {d.A, = E E §'^(^^(0, ^.(0) (5) 
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It must be noted that for notational convenience, we drop the time dependence of the functions 
I{N, {dij},l), IcXN, {dij},l) and Ic^^kiN, {dij},l) following the convergence of the algorithm 
or whenever the spatial configuration is fixed over time, and denote them by I{N, {dij}), 
Ici{N, {dij}) and /c.,fe(A^, {dij}), respectively. We also denote the aggregate interference of the 
worst-case scenario, optimal scenario and that of the output of the algorithm by Iw{N, {dij}), 
Io{N, {dij}) and laiN, {dij}), respectively. We can now define the main algorithm: 

Main Algorithm: Clusters scan all the frequency bands 61 , ■ ■ ■ , 6^ m an asynchronous manner 
over time. Each cluster chooses the frequency band in which it experiences the least aggregate 
interference from other clusters. In other words, at time I, a cluster, say Ci, updates its state 
according to the following rule 



where Sj(/ + 1) is the new state of Ci updated at time I. 

For this purpose, the cluster head scans all the frequency bands and estimates/measures the 
interference it experiences in each frequency band. The cluster head chooses the new transmission 
frequency band according to Eq. the decision criterion in the Main Algorithm. 



A. Convergence 

Theorem 4.1: Given any reciprocal channel model, the Main Algorithm converges to a local 
minimum in polynomial time in A^. 



B. Performance Bounds 

Theorem 4.2 (Upper Bound): Let Ia{N, {dij}) denote the aggregate interference of all the 
clusters corresponding to the state of the algorithm following convergence (see Theorem 14.11) . 
and IwiN, {dij}) to be the aggregate interference for the worst case interference scenario (all 
clusters transmitting in one frequency band), then 



Si{l + 1) = argmin/c,,fc(A^, {dij}, I) 



(6) 



IV. Main Results 



Proof: The proof is given in Appendix HI 



□ 



Ia{N,{di,})<lUN,{d,,}) 



Proof: The proof is given in Appendix III-A[ 



□ 
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In order to obtain further performance bounds for the Main Algorithm, one must analyze its 
behavior following convergence, over a configuration induced by the spatial configuration of the 
network nodes. However, such an analysis for the general case is non-trivial. We therefore focus 
our attention on a specific class of spatial configurations, where all the clusters are co-linear (lie 
on a line). This is referred to as a Linear Array. For a given A^, we assume that the clusters are 
located in [0, (A^ — l)d], where d is a constant. 

Definition 4.3: A Uniform Linear Array of clusters is a Linear Array in which (ij j+fc = kd, 
for 1 <i,i + k < N. 

Theorem 4.4 (Optimal Strategy): Let Io{N, {dij}) be the aggregate interference of the optimal 
strategy for a given linear array of clusters located in [0, (A^ — l)d]. Then, 

\im ^To(N,{diA)> —2C(v)^ (7) 
N^oo N ^ ^ ~ r'^ ^ d-^ 

where ({r]) is the Riemann zeta function. 

Proof: The proof is given in Appendix IH-BI □ 

Theorem 4.5: Uniform Linear Array achieves the bound in the statement of Theorem 14.41 as 

A^ oo. 

Proof: The proof is given in Appendix IH-CI □ 
Theorem 4.6: If r = 2 and > 2 in the statement of Theorem 14.51 then the optimal strategy 
is the alternating assignment for any A^. 

Proof: The proof is given in Appendix III-D[ □ 
We can use Theorems 14.41 and 14.51 to upper bound the performance of the algorithm as given 
by the following Corollaries: 

Corollary 4.7: Let /a(A^, {dij}) denote the aggregate interference corresponding to the output 
of the Main Algorithm and /o(A^, {dij}) be the aggregate interference of the optimal strategy. 
Then, for any spatial configuration of clusters in [0, (A^ — l)d\ we have 

as A^ ^ OO. 

Proof: The proof is given in Appendix III-EI □ 
Corollary 4.8: If the array in the statement of Corollary 14.71 is a Uniform Linear Array, then 
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(9) 



UN.d) ^ 
UN,d) - 

where Ia{N,d) = Ia{N, {dij}) and Io{N,d) = Io{N,{dij}) for a Uniform Linear Array in 
[0,(iV-l)rf]. 

Proof: The statement of the Corollary follows by letting d^ax = f^mm = in Corollary 

m □ 

C. Time-varying Setup 

For cluster q, i = 1, 2, ■ ■ ■ , A^, we consider an activity indicator state ai{l), such that aj(/) = 1 
and ai(l) = correspond to being active and inactive at time /, respectively. Let Pq{1) and 
Pi' {I) be the probability of Cj being in activity indicator state and 1 at time /, respectively. 
The evolution of the probabilities is given by: 

Po^'(/ + l) a Pod) 

Pf(/ + 1) j \l~a a ) [p^'il] 

which corresponds to a symmetric two-state Markov model. We assume P^'^O) = 1 for all q, 
z = l,--- ,iV. 

In this case, the convergence of the algorithm is not guaranteed. Furthermore, the upper bound 
on the performance of the algorithm does not hold in presence of time-varying statistics. In order 
to further analyze the performance of the algorithm, we use the following additional assumptions: 

• The clusters update their frequency band asynchronously according to the same temporal 
statistics. 

• The update process is modeled by a Poisson process of rate i.e., each cluster updates 
its frequency band with a rate j^^- 

• The number of users switching on/off in each time slot is much less than the total number 
of users (near equilibrium scenario). 

D. Dynamics of The Algorithm 

Let /(t) denote a continuous-time approximation to the aggregate interference of the network 
at time t. It can be shown that 
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« = -tiEm] - m) (10) 

where 8 denotes the ensemble average (over different update patterns), r = NAT and p is 
geometrical constant showing the effective number of interacting neighbors to a cluster including 
itself (the details are given in Appendix IIII-A|) . 

Eq. (fTOl) is an approximate differential equation describing the behavior of the algorithm near 
equilibrium. In other words, the aggregate interference decreases exponentially with rate p/r 
near the equilibrium point. 

We can model the change in the number of active clusters by two Poisson counters of rate 
A [5]. We associate two independent Poisson counters dN^ and dN_ to the users going on and 
off, respectively. We have 

E{dN+) = E{dN_) = Xdt (11) 

Each cluster when activated, approximately experiences the instantaneous normalized aggre- 
gate interference of the network, due to spatial ergodicity. If we define I(t) = £[I{t)] and 
Xa(t) = £[Ia{t)], under the assumption of A being small compared to ^, we have the new 
dynamics in the Ito form [7] as 

dl{t) = --(J(t) - Ia{t))dt + ^I{t){dN+ -dN_) (12) 
r N 

Under this model, in the steady state the variance settles down to 



{m - E[x{t)]y 



given 

(please see Appendix IIII-BI for details). 

According to our model for cluster activities give by Eq. Q, A = A^^(l — q)/2t. Thus, we 
get the following trade-off inequality 

8(1 - a) 



P 

in order to have a finite variance in the steady state. 



< 1 (15) 
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E. Discussion of the Results 

Theorem 14.11 guarantees the convergence of the Main Algorithm in polynomial time in A^, 
regardless of the configuration of the nodes in the network. The first performance result is stated 
in Theorem 14. 2[ which gives an upper bound on the performance of the algorithm. This result 
is independent of the topology of the network and is a direct consequence of the structure of 
the Main Algorithm. 

In order to obtain further performance bounds for the Main Algorithm, as stated earlier, we 
have chosen the class of Linear Arrays for our analysis. Theorem 14.41 is the main result on the 
asymptotic performance of the Main Algorithm for Linear Arrays. The Riemann zeta function in 
the statement of Theorem is merely a consequence of the path loss model for the transmission 
and the fact that at each time slot, only one transmitter and one receiver are active in each 
cluster. However, the achievability of the bound in Theorem 14.41 is not trivial, since the optimal 
strategy of frequency allocation is not known for a general spatial configuration of the clusters. 
Theorem 14.51 proves the asymptotic achievability of the bound in the statement of Theorem 14.41 
for Uniform Linear Arrays. Theorem 14. 5[ states that the alternating frequency band allocation 
achieves the bound in Theorem 14.41 as ^ oo, and is therefore the optimal strategy of frequency 
allocation. In order to analyze the optimal strategy for finite N , we focus our attention on the 
case of r = 2, i.e., two frequency bands. Theorem 14. 6 [ states that the optimal strategy is the 
alternating frequency allocation for finite A^, when r = 2. 

Corollaries 14.71 and 14.81 compare the performance of the Main Algorithm with respectively that 
of the optimal strategy for Linear Arrays and Uniform Linear Arrays. This is done by combining 
the upper bound of Theorem 14.21 with the lower bound of Theorem 14.41 The bound in Corollary 
14.71 is not as tight as the original bounds, since it is worst-case and applicable to any spatial 
configuration of clusters. As the simulation results show in next section, the algorithm performs 
significantly better than these bounds. 

The stochastic analysis for the time-varying case results in a simple trade-off inequality for 
design purposes. In other words, for a given time-varying statistics, we can design the update 
rate to guarantee a finite steady state variance. 

The geometrical parameter p can be empirically estimated for different network topologies. 
However, theoretical estimates are possible. For instance, for a uniform linear array of clusters. 
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each cluster has two effective neighbors. Therefore, an estimate of p = 3 seems reasonable and 
is verified by simulation results (Please see Section IV]). If A = 0{N^^'') for some e > 0, as 
N ^ CO, the inequality (fT4)) always holds. Therefore, the algorithm converges in both mean and 
variance in the sub-linear regime. 

F. Extension of the Results 

Generalization of the result of Theorem 14.61 to r > 2 is not straightforward, since the 
combinatorial possibilities of the assignments grow exponentially with r. Furthermore, gen- 
eralization of all the above results to higher dimensions is a non-trivial problem. The reason is 
that in dimensions greater than 1, the degrees of freedom for the cluster interactions increase 
dramatically. We are currently studying this scenario extensively. 

As noted earlier, any model with channel reciprocity suffices for the convergence of the Main 
Algorithm (See the proof in Appendix U). Moreover, the upper bound on the performance of 
the algorithm given in Theorem 14.21 holds for any model with channel reciprocity. The other 
performance bounds, in their current format, rely on the specific path loss model. However, it 
may be possible to generalize the same method of analysis to other reciprocal channel models. 

The stochastic analysis can be generalized to other statistical models for the activity of the 
clusters over time. The generalization of the results to r > 2 is straightforward. 

V. Simulation Results 

Fig. 2 shows the performance of the algorithm on different arrays of 100 clusters in one and 
two dimensions. Fig. 2 (a), (b) and (c) show the normalized aggregate Shannon capacity (both 
the optimum/near-optimum value and that of the output of the algorithm) and the normalized 
aggregate interference of the network as a function of time for a uniform linear array (with r = 2), 
integer Z2 (rectangular) lattice (with r = 4) and a hexagonal lattice (with r = 4), respectively 
(Please see [3]). The normalized aggregate Shannon capacity is defined as the aggregate Shannon 
capacity divided by the number of clusters. Similarly, the normalized aggregate interference is 
defined as the aggregate interference divided by the number of clusters. Here d = I, Pq = I and 
1] = 2. For the rectangular and hexagonal lattices, the computation of the optimum frequency band 
assignment is very complicated and finding it by exhaustive search was beyond the capabilities 
of our simulation platforms. Instead, we have compared the performance of the algorithm to that 
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of the 1 : 4 frequency reuse pattern as a near-optimal candidate. As it can be observed from the 
figure, the minimization of the aggregate interference results in an overall increasing behavior of 
the capacity. In all cases more than 90% of the capacity of the optimal (near-optimal) centralized 
frequency assignment is achieved. 

Fig. 3 shows the performance of the algorithm on a uniform linear array along with the lower 
and upper bounds we have obtained in Section IIV-BI Here r = 2, d = 1, Pq = 1 and rj = 2. 
For the initial condition of the algorithm, we let all the clusters to be in frequency band hi. 
The updates are repeated until the convergence is achieved. As we observe from the figure, the 
algorithm performs significantly better than the upper bound we have obtained and is less than 
1 dB away from the alternating assignment and the lower bound. 

Fig. 4 (a) and (b) show the performance of the Main Algorithm, the worst case and the 
1 : 4 frequency reuse pattern for rectangular and hexagonal arrays of clusters, respectively. Here 
d = 1, r = A, Pq = \ and r] = 2. Although we have not established any lower performance 
bounds for two-dimensional arrays, the algorithm performs very closely to the centralized 1:4 
frequency reuse solution. 

In Fig. 5, the dynamics of the algorithm for a uniform linear array of 100 clusters is shown. 
Here a = 1 which corresponds to the case of no time-varying statistics. The empirical curve is 
averaged over 500 different update patterns. The theoretical curve corresponds to an exponential 
with rate —p/r for p = 3 and r = 1. As it can be observed from the figure, the theoretical 
estimate matches the empirical curve very well. 

Fig. 6 shows the theoretical and empirical steady state variance of the aggregate interference vs. 
switching rate (1 — a) for a uniform linear array of 100 clusters. The empirical curve is obtained 
by averaging over 500 different realizations of the update process. As it can be observed from 
the figure, the theoretical variance (with p = 3) matches the empirical variance. 

VI. Conclusion 

We have proposed a distributed algorithm for finding a sub-optimal frequency band allocation 
to the clusters in a network. We have also derived some performance bounds for the special 
case of linear arrays of clusters. Simulations prove that the algorithm performs significantly 
better than the performance bounds we have established. We have also derived a stochastic 
differential equation describing the behavior of the algorithm in presence of time-variation in 
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the activity of clusters near the equilibrium point. A trade-off inequality to guarantee stability in 
the performance of the algorithm is established. The stochastic modeling framework opens the 
possibilities of both open loop and closed loop stochastic control. These problems are currently 
being studied. 

Appendix I 
Proof Of the Convergence Result 

First we prove the following Lemma: 

Lemma 1.1: I{N, {dij}, I) is a non-increasing function of / that is bounded from below. 
Proof: First of all, we need to show that I{N, {dij}, I) has a lower bound. We know that the 
aggregate interference of all the clusters is a non-negative quantity. Therefore, I{N, {dy}, I) >0 
for all /. 

Secondly, we need to show that I{N, {dij}, I) is a non-increasing function of /. Without loss 
of generality, we assume that q has been transmitting in band bj. I{N, {dij}, I) can be written 
as 

I{N, {d,,}, l) = J2 E ^-^r.,k{N, {d^,}, I) + 2h^,{N, {rf,,}, /) (16) 

where we have used the channel reciprocity. After the update, the algorithm implies that q 
chooses the new band hk according to the decision criterion Q in the statement of the Main 
Theorem. Therefore, 

I{N,{d,,},l + l) = J2 Yl Ic^AN,{d^J}J) + 2h^,k{N,{d,,},l) (17) 
= I{N, {rf,,}, /) + 2(4,,(iV, {d,,}, I) - 4,,(iV, {d,^}, I)) 

where k satisfies the decision criterion According to Eq. we have Ic^^k{N, {dij}, I) < 
Ic-j{N, {dij}, I), for all j ^ k. Therefore, 

I{N, {di,}, l + l)- I{N, {d,,}, I) = 2(4,fc(iV, {d,,}, I) - 4,,(iV, {d,,}, I)) < (18) 

which gives the statement of the Lemma. It must be noted that channel reciprocity is sufficient 
for this proof to hold. □ 
Proof of Theorem \4.1\ 

According to Lemma [LTI I{N, {dij}, I) is a lower bounded non-increasing function of /. It can 
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take at most r distinct values, corresponding to the r different frequency band assignments 
to the clusters. Therefore, 3/o such that V/ > Iq we have I{N, {dij}, / + 1) = I{N, {dij}, I) and 
no cluster updates its state up to an isomorphism from {0, 1, ■ ■ ■ , r — 1} to itself. 

After each update, the change in I{N, {dij}, I) is at least of order — y^, where dmax = 
maxij{dij}. The aggregate interference of all the clusters is of order A^. Therefore, we need 
at most 0(A^''+^) switches to reach the final configuration of the algorithm. After any round 
of updates, during which all the clusters have updated their state, at least one cluster changes 
its frequency band. Therefore, the total number of updates is 0{N^^'^), which is bounded by a 
polynomial in A^. Therefore, the Main Algorithm converges to a local minimum in polynomial 
time in A^. 

Appendix II 
Performance Bounds Proofs 

A. Proof of Theorem \4.2\ 

Lemma 2.1: Let Ii{N, {dij}) = w '-^^ interference on q when all the other clusters 
are co-band with it. Then, 

h^,k{N,{d,,},l) < }M{N,{dij}) 
Proof: At time /, q chooses a frequency band, say k E {1, 2, ■ ■ ■ , r} such that Ic^^kiN, {dij}, I) < 
Icij{N, {dij}, I), for all j ^ k. Therefore, 

rh^AN, {d^,}, 1)<Y1 ^-«.^(^' = ^^(^' {d^J}) (19) 

which gives the statement of the Lemma. □ 
Proof of Theorem \4.2\ 

Let IcXN-,{dij}) denote the value of /c.^fc(A^, {dj^}, /) following convergence. Using the result 
of Lemma 12. 1[ the aggregate interference can be written as 

h{N, {d,,}) = J2 IcAN, {di,}) <^Y1 ^^^^-^ = {^^^•}) (20) 

i i 

which gives the statement of the Theorem. 
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B. Proof of Theorem \4.4\ 

Lemma 2.2: Let Iw(N,d) denote the aggregate interference of the worst case scenario for a 
uniform linear array of N clusters in [0, (N — l)d]. As — > oo, we have 

l-I^{N,d) ^2av)^ (21) 

where (^(■) is the Riemann zeta function. 
Proof: We have 



1 1 P ^ ^ 1 

*=i i=i ' ' 

N oo 



(22) 



< i^y2yl = 2C(r/)^ 

N d'^ ^ ^ j-n " rf'? 

i=i i=i 

Let Ii(N,d) = /j(A^, {djj}) denote the interference experienced by the cluster q located at 
(i — l)d, when all the other clusters are co-band with it, for a uniform linear array of N clusters 
in [0, (A^ - l)d]. We can write 

We need to show that for any small enough e > 0, 3A^o such that for A^ > Aq^ 2^(77)^ — 
^/^(AT, d) < ef. Let's take N such that JnT^ < f - For N > 2{N + 1) we have 

-I^{N,d) = -J2l,{N,d)>- J2 ^^(^'^) (24) 

«=1 j=JV+l 

Clearly for all iV < i < A^ - N, 

2C(r^)^ - /.(AT, rf) < 2^ y -J— < (25) 

Therefore, 

^/»('V.<') > '^<('V,<') > ^§EST'(2C(.,) - f) = ^i«(2C(.) - i) (26) 

i=N+l 

If we choose A'o large enough so that A({r])^^ < e/2, for all A^ > A'o we will have 

i^^^^|.4CM^ < ./2-../2^. (27, 
which proves the statement of the Lemma. □ 
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Lemma 2.3: For a given distribution of clusters in [0, {N — l)d], we have 

UN,{d,,}) ^ 

as — > oo. 

Proof: Let's assume that N is even. Suppose that we fix the first cluster at and the 
Nth cluster at {N — l)d. We define a vector interaction field between any two clusters as 
Fij — a/dj^^ with a — Pq/t], and suppose that each cluster is denoted by a point on the interval 
[0, (A^ — l)d] interacting with the others according to our vector field. Therefore, if we let the 
points be initially distributed on the line with distances {rfy }, the final equilibrium configuration 
will be when Iyj{N,{dij}) is minimized. This is because Iw{-,-) acts as a potential function 
for the system due to the definition of the interaction vector field. Let {o?y (/)} denote the final 
distances of the clusters following equilibrium. If we consider /^(-j •) as a function of adjacent 
distances Aj = dj^j+i for i = 1, • • • ,N — 1, it will have the following form: 

W{A.}) = 2 5: (28) 

i,j:j>i \2^i V 

which is clearly a convex function of A^, i = 1, 2, • • • , — 1. To find Aj for i = 1, 2, ■ ■ ■ , N — 1 
in the equilibrium configuration, we need to minimize Iyj{N,{Ai}) subject to the constraint 
Aj = Nd. Since, the constraint is linear in Aj for i = 1, 2, • • • , — 1 and the second 
derivative of Iiu{N,-) with respect to any Aj, i = l,2,---,A'^ — lisa sum of some positive 
terms, there is no i, for which the second derivative of Iyj{-, •), the objective function, is zero. 
Therefore, the minimization has only one unique solution which is the equilibrium state (This is 
intuitively clear from the physical model of a number of points interacting on a one dimensional 
lattice). Since the system is in equilibrium under the distribution given by {c?ij(/)}, we have 

D(l/(Ar,{A,}) 

where D{-) denotes the exact differential and Afc(/) = dk,k+i{f)- We claim that if we choose 
Ak — d for all k, then for every e small enough, we can find an TV so that for all N > N we 
have: 

This means that we can get arbitrarily close to the equilibrium point of the system. Since, the 
objective function is continuous and differentiable with respect to Aj, for i = 1, 2, • • • ,N — 1, 



= (29) 

Afc=Afc(/) 
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we have chosen to approach the equilibrium point uniformly with respect to all Aj, for i = 
1,2, ■ ■ ■ , N — 1, i.e., choosing the same amount of variation for all the coordinates. From the 
properties of a potential function for a vector field, we know that 

V d^N, {A J) -1vf<^Vifi nn 

^ dK, A.^A.(/) -nI^^^'^n1^\^^\ (31) 

where Fi is the aggregate interaction field of all the other particles on the zth particle. If we set 
Ai = d for all the users, then 

N-i+l i 

\F'\ = -Ejjsr^<^{ii^^^)-Ej^) (32) 

j=i+i ^-J I j=\ J 

for all % < N/2, and = |-F/v-,;+i| for i > N/2. Therefore, we can bound the expression in 
the righthand side of Eq. (|3T1) as 

1 2 ^^"^ ' 1 2 ^'"^ 

N E l^^l < E + 1) - E -;^) = ^ (]V E (33) 

i i=l j = l i=l 

where Oj = C('7 + 1) ~ X]j=i jh^- Clearly, Ymvi^ooCa = 0. Therefore, if we choose N so that 
Yl^N+i < similar to the proof of Lemma 12.21 we can choose large enough to 
make sure that Eq. (|30l) holds. This means that the uniform spatial configuration asymptotically 
coincides with the global minimum of the system, which is the optimal spatial distribution of 
clusters with least aggregate interference. Since the result holds for large A^, fixing the first and 
last cluster at and {N — l)d does not affect the result. The case for odd can be treated 
similarly, since excluding one of the users does not change the asymptotic result. □ 
Lemma 12.31 gives an important result: the worst-case aggregate interference of any spatial 
distribution of N users in [0, (A^ — is greater than or equal to that of a uniform linear 
array of A^ clusters in [0, (A^ — l)d], for sufficiently large A^. 

Lemma 2.4: Let Ck{N, {dij}) denote the set of clusters in frequency band hk corresponding to 
the optimal assignment strategy for a given spatial configuration of A^ clusters in [0, {N—l)d\. Let 
Nk{N,{dij}) denote \Ck{N,{dij})\. Then, \^Nk{N , {dij})^ ^^^^ is an unbounded sequence, 
for all /c = 0, ■ ■ ■ , r — 1 and any spatial configuration of A^ clusters in [0, (A^ — l)d]. 

Proof: Let5(Ar, {rf,,}) = {Nk{N, {d,j})\k = 0, ■ ■ ■ ,r-l] ^ndN*{N, ^ iniS{N, {d.^}). 

Let's assume that {N*{N,{dij})}N=i^2,--- is a bounded sequence. That is, there exist numbers 
N and M and a specific spatial configuration of the clusters for any N > N, such that 
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N*{N, {dij}) < M. Let N*{N, {dij}) correspond to frequency band bk* (whereas dependence 
on is implicit). We have 

iv.w,}) = lE E ^ + ^ E ^ (34) 

k^k* c,,c,eCk{N,{d,,}) *J c„Cj&Ct,*iN,{d,j}) *J 

where we use the additional argument of r to show the implicit dependence of the aggregate 
interference on the number of frequency bands. We take the same ensemble Cj for i = 1, ■ ■ ■ , 
and assign any configuration of frequency bands bk for k ^ k* to all q G Ck*{N, {dij}) and keep 
the frequency assignment of the rest of the clusters as in the optimal strategy. Let s,; denote the 
index corresponding to the new frequency band assigned to q G Ck*{N, {dij}). The normalized 
aggregate interference of this new ensemble will be 

l/(r-l,Ar,{d.,}) = 1^ Yl ^ (35) 

k=ik* C„C,€Ct:(N,{d,j}) 



N ^ ^ ^ d: 

k^k* c,gCfc*(Af,K,}) CjGCfc(7V,K,}) 



We also know that /(r - 1, A^, {dij}) > Io{r - 1, A^, {dij}), where /^(r - 1, A^, {dij}) is the 
aggregate interference of the optimal frequency band assignment to the same ensemble, using 
r — 1 frequency bands. We define 

Si^^id,})^ ' J: ^ (36) 

m{d,})^ ' E E ^^(^^'^) (37) 

Then, we have 

^Io{r,N,{d,,}) = l/(r - 1, N, {d,,}) + (g(N, {d,,}) - h{N, {ci.,}))(38) 

> ^/o(r - 1, AT, {d.,}) + ^^^^^^^{g{N, {rf,,}) - hiN, {d,,})) 

Both g{N, {dij}) and /i(A^, {(iij}) are upper bounded by 2({r])^, where 5 is the minimum 
distance between two clusters, based on our assumptions. We have ^ ^^j^'^'^^^ \ {g{N^ {dij}) — 
h{N, {dij}))\ < A({ri)^[^). Therefore, we can choose A^ sufficiently large with a specific 
corresponding ensemble, in such a way that ^ \g{N, {d^j}) — h{N, {dij}) \ becomes arbi- 

trarily small. Equation (l38l) states that the normalized aggregate interferences using r frequency 
bands and r — 1 frequency bands for a specific configuration may become arbitrarily close. 
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This is clearly not possible. Suppose that we have r — 1 frequency bands bQ,bi, ■ ■ ■ , &r-2- 
Consider the optimal frequency band assignment to the clusters q for i = 1, 2, ■ ■ ■ , iV in [0, (A^ — 
l)d]. Since ^k^k{N, {dij}) = N, there exists aO<ko<r-2 such that Nka{N, {dij}) > 
N/ (r — 1). If we consider the clusters in Cko{N, {dij}) and allow the assignment of an additional 



frequency band 6^-1 to them, then according to Theorem 14.21 there exists a frequency band 
assignment using the two bands bk^ and br, to the clusters in Cko{N, {dij}), for which the 
normalized aggregate interference of all clusters in Cko{N, {dij}) is at least half of that in the 
worst case scenario. Therefore, adding a frequency band br-i, decreases the normalized aggregate 
interference by at least 

2iV 5Z ^ ^ (r-l)^+^^^'^^dv ^^^^ 

for sufficiently large A^, according to Lemma [23l Therefore, j^Io{r — ^,N,{dij}) can not be 
arbitrarily close to j^Io{r,N, {dij}) for any spatial configuration of clusters, as ^ oo. Thus, 
{N* {N, {dij})} N=i^2,-- is an unbounded sequence. Since N*(N, {dij}) = ini S{N, {dij}) and 
S{N, {dij}) is a finite set, we conclude that {Nk{N, {(ijj})}Ar=i,2, ■■ is an unbounded sequence 
for all A; = 0,1,- ■■ ,r - 1. □ 
Lemma 2.5: Let Nf^ denote \Ck\ corresponding to the optimal strategy for a given spatial 
configuration of clusters in [0, (A^ — l)d]. Then, we have ^'^^^'^''^''''^^^^^'•''j^'^k ^ k, for 

sufficiently large A^. 



Proof: Using the result of Lemma 12.31 we can write 

1 1 A^ 

—urn, {rf.,})|c„c,ec, > -i^UNk, —d) (40) 

since the clusters are located in [0, (A^ — l)d] and A^^ is unbounded for all k (Lemma 12.41 ) as 

A^ ^ oo. 



□ 



Proof of Theorem \4.4\ We have 



r— 1 

1 1 



lim -/,(Ar,K,}) = Jiin^^E E UNk,{dim}) (41) 



Af^oo A^ 



1 '■"-^ v^r-i prr, 

k=0 
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where we have used Lemmas |2.2[ 12.41 and l2.5[ The expression Xlfe=o ^fc^^ achieves its minimum 
when Nk = N/r, for all A; = 0, ■ ■ ■ , r — 1. Therefore, 

Ita -UN, {d,}) > ^-%^UN, d) > - Jta -UN, d) = -2C(r,)^ (42) 
which gives the statement of Theorem I4.4[ 

C. Proof of Theorem 14.51 

Lemma 2.6: Let Io{N,d) = /o(A^, {c^ij}) and Iw{N,d) = /u;(A^, {c^ij}) denote the aggregate 
interference for the optimal alternating strategy and worst-case scenario, for a uniform linear 
array in [0, (A^ - l)d]. Then, 

rI^{[N/r\ , rd) < UN, d) < rI^{\N/r'\ , rd) 
Proof: We have a uniform linear array of clusters in [0, (A^ — In case of the 
alternating frequency assignment, the clusters corresponding to the hk band make a uniform 
linear array of [A^/rJ < Nk < \N/r~\ clusters each distance rd apart from its neighbor(s), for 
all k = 0, ■ ■ ■ ,r — 1. Therefore, the aggregate interference is equivalent to the sum of aggregate 
interferences of the r disjoint sets of users corresponding to the r frequency bands. Clearly, all 
the clusters in each subset have the same frequency band and therefore, the aggregate interference 
for each subset will be Iw{Nk, rd). We can write Io{N, d) as Io{N, d) = Y7k=\ ^wi^k, rd). Since 
Iw{N, ■) is an increasing function of A^, we can sandwich Iw{Nk,rd) as 

I^{YN/r\,rd) < lUNk,rd) < lU\N/r],rd) (43) 

Thus, 

r-l 

rlU[N/r\ ,rd)<J2 UNk, rd) < rU\N/r] , rd) (44) 

k=0 

which gives the statement of the Lemma. □ 
Proof of Theorem 14.51 

Using Lemmas 12.61 and 12. 2[ and using the Sandwich Theorem in Calculus we can see that, 

lim ^Io{N,d) = lim l-rI^{\N/r],rd) = lim M!l2C(r/)-^ = -2C{r])^ (45) 
from which the statement of the Theorem follows. 
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D. Proof of Theorem \4.6\ 

Lemma 2.7: Let's consider a specific assignment of the frequency bands to clusters in a 
Uniform Linear Array, such that there are 3 successive co-band clusters. Switching the frequency 
band of the middle cluster decreases the aggregate interference, regardless of the frequency bands 
of the rest of the clusters, if C{ri) < 2, where C(-) is the Riemann zeta function. 

Proof: Without loss of generality, assume that Cj_i, c, and q+i are in frequency band hi. 
The interference of all the other clusters on Cj is clearly greater than or equal to 2PQ/d^. If it 
switches to band bo, the interference of all the other clusters on q will be less than or equal to 
2(C(^) ~ ^)Po/d^, for which the equality occurs when there are infinite number of users other 
than Cj„i and Cj+i in band Bq. Therefore, the aggregate interference will decrease after switching 
if 2(^(77) — 1) < 2, which gives the statement of the Lemma. □ 

Lemma 2.8: For any given frequency assignment, we can find a configuration which has 
smaller aggregate interference, and can be partitioned into subsets of clusters with the property 
that in each subset all the clusters have alternating frequency assignment. 

Proof: Using Lemma I2.7[ we can break any 3 or higher successive co-band clusters into 
at most 2 successive co-band clusters, by repeatedly switching the frequency band of the middle 
cluster. We will then get a configuration which has smaller aggregate interference compared to 
the original one (according to Lemma 12.71) . We can therefore partition this new configuration 
into subsets of the stated property, by just putting the boundaries midway between every two 
successive co-band clusters, i.e., if we have k pairs of co-band clusters, we can partition the 
clusters into k + 1 subsets 5'i,5'2, ■■■ ,5*^+1, in such a way that all the clusters in Si have 
alternating frequency assignment, for alH = l,2,---,A; + l. □ 

Lemma 2.9: Let's consider the subsets S*!, 5*2, ■ ■ ■ , Sk+i given by the partitioning procedure in 
Lemma [Z8l Then, switching all the frequency bands in Si, decreased the aggregate interference. 

Proof: Without loss of generality, let's assume that there are Ns-^ = 21 clusters in Si and 
the first cluster is in band bo. Since A^Si is even, the last cluster in Si is in band hi. Also, 
the first cluster in 5*2 IJ 'S'3 U ' ' ' U ^k+i is in band 61. Let Ii{b) and /i(a) respectively denote 
the interference of all the clusters in C\Si on the clusters in Si in band 61, before and after 
switching. Therefore, 

Ah 4 Ii{b) - Ilia) >5(l-^ + :^-i + -- - + 721V - (2IF) (46) 
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Similarly, let Io{b) and Io{a) respectively denote the interference of all the clusters in C\Si 
on the clusters in Si in band bo, before and after switching. Therefore, 

|A/o|^|/o(&)-/i(a)| < ^{(C(r])-l)-(C(^)-l-2^) + (C(^)-l-^-^) (47) 

2v 317 4»7 ) + (CW-l-jj lif;)} 

dv\2v^4'^ (2Z)V 
Therefore, the difference in the aggregate interference of the system before and after switching 
will be 

A/ 4 2(A/, + A/„) > 2^{(1 + + + + i + + (48) 

> 2§{1-1C(,)} 

1 — ^Civ) is clearly positive for all i] > 2. This means that switching all the frequency bands 
in 5*1 decreases the aggregate interference, i] > 2 also satisfies the bound in Lemma [277] Thus, 
the optimal solution is the alternating solution for > 2. The case for odd N can be treated in 
a similar manner. □ 

Proof of Theorem \4.6\ 

Suppose that we are given a specific frequency band assignment for N clusters in a Uniform 
Linear Array. Using Lemma [Z8l we can partition the clusters into subsets ^i, ^2, ■ ■ ■ , S^+i- Then, 
using Lemma |239l we switch all the frequency bands in Si. After switching, the last cluster in 
5*1 will not be co-band with the first cluster in S2- Therefore, we can combine Si and 5*2 to 
get 5*2 = 5*1 IJ 5*2, reducing the number of partitions to k. We can now switch all the frequency 
bands in ^2 and combine it with S^ to get ^3 = S[\JS2 = "S*! |J 6*2 |J 6*3. If we repeat this 
procedure a total of k times, the resulting configuration is clearly an alternating assignment. Since 
the aggregate interference decreases in each of the k steps, the final alternating configuration 
has smaller interference than the initial configuration. Noting that the initial configuration was 
arbitrary, we conclude that the alternating frequency assignment strategy is the optimal strategy 
of frequency assignment for a finite Uniform Linear Array with r = 2 and 77 > 2. 

E. Proof of Corollary \4. 71 

Clearly, we have Iw{N,dmax) < Iw{N,{dij}) < In,{N ,dmin), where dmin and dmin are the 
minimum and maximum spacing between the clusters, respectively. The statement of Theorem 
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Io{N,{dij}) ^ ]_/drniny ^^^^ 



I4.4[ can be written as 

Iw{N, drain) ~ r"? V d 

as iV ^ oo. We also have Io{N,{dij}) > Io{N,dmax), where Io{N,dmax) is the aggregate 
interference of the optimal frequency band allocation for a Uniform Linear Array in [0, (A^ — 



^)dmax]- Using Theorem 14.51 we have 



Af— >oo I^(^N, drain 

Combining Eqs. (|49] ) and (l50l) we get 



1. -^oi-^j dmax) 1 / dnaa \ /rrw 

hm = — (50) 



max 



Io{N, {dij}) ^ 1 / dra.a , ^^^^ 



Iu,iN,dmin) ~ r'^\mm{dn,ax,d} 
Using Theorem Sill a bound of IaiN,{dij}) < ^I^uiN, {dij}) < ^I^iN^dmin) is obvious, 
where Ia{N, {dij}) is the aggregate interference corresponding to the output of the Main Algo- 
rithm. Combining this bound with Eq. (|5TI) . we get 

as N oo. The latter gives the statement of the Corollary. 

Appendix III 
Dynamics Of The Interference 

A. No Time Variation In The Activity of Clusters 

We consider the case of two accessible frequency bands, and hi, for mathematical con- 
venience. Generalization to more than two frequency bands is straightforward. We associate 
ej = — 1 and ej = 1 to clusters in band 6o and hi, respectively. 

If Cj is active at time t, it will experience an interference 

Ic.{t) = \{h + emY.'-^) (53) 

where Jj is the worst case interference experienced by cluster q. 

We define the band hj to be appropriate for cluster q, if q is assigned in band hj in the 
target configuration of the algorithm. If Cj is not assigned in its appropriate band at time t, it 
will increase the aggregate interference by 
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2A4(t) = 2|5^!^ 



(54) 



If a total of M(t) users Cj, z G {ii, ■ ■ ■ , «M(t)} are not in the appropriate frequency bands, the 
aggregate interference will be 

M{t) 

(55) 



k=l ■ j^ik 

where /„ is the target performance of the algorithm. 

Assuming spatial ergodicity, we can average over an ensemble of different update patterns. 
Thus, 



6,(t) 



S[I{t)]-8[h]=2M{t)8\^Y. d"^. 

where 8 denotes the ensemble average. 

After any update, the average change in 8[I{t)], will be 



(56) 



A.[/(t)]^^24i$: 



(57) 



where p is a geometrical constant showing the effective number of interacting neighbors to a 
cluster including itself. We are assuming that every cluster which is not in the appropriate band, 
increases the interference on an effective number of its neighbors (this is a linearization near 
the equilibrium point). 

Combining Eqs. (l56l) and (l57l) . we get 

^ = -^M/Wl-£|/J) (58) 

where AT is the inverse update rate. On the time scale of the updates, using the ansatz r = NAT 
one can write 

« = -t(sm\ - (59) 
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B. The Variance Equation Under Time-varying Activity of Clusters 
Using the Ito rule, the variance equation associated with Eq. (fT2)) is 

dl\t) = -'^I{t){l{t) - la{t))dt (60) 

r 

+ ((X(t) + ^m)" -l\t)yN^ + {(X{t) - |x(t))' - J^(t))rfiV_ 
By simplifying the expression and taking expectations we get 

Therefore, in the steady state the variance settles down to 



given 



{i{t)-Em]y 



j-2 '^PI'T t2-t2 ^^^^ (-62^ 

"2p/r-32A/iV2 " "iV2p-16Ar ^ ^ 
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Fig. 2. Normalized aggregate Shannon capacity and normalized aggregate interference curves for arrays of 100 clusters vs. time, (a) uniform 
linear array, (b) rectangular lattice, and (c) hexagonal lattice. 
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Fig. 3. Normalized aggregate Shannon capacity and normalized aggregate interference curves for uniform linear arrays vs. N 




Fig. 4. Normalized aggregate Sharmon capacity and normalized aggregate interference performance curves for 2 dimensional arrays vs. N. 
(with r = 4) (a) rectangular lattice, and (b) hexagonal lattice 
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Fig. 5. Normalized aggregate interference vs. time for a uniform linear array of 100 clusters with a = 1. 



0.018 




0.002 1 ' ' ' ' ' ' ' ' 1 

0.001 0.002 0.003 0.004 0.005 0.006 0.007 0.008 0.009 0.01 
Switching Probability (1 -a) 



Fig. 6. Normalized steady state variance vs. switching rate for a uniform linear array of 100 clusters. 
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